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Many experiments in quantum information aim at creating multi-partite entangled states. Quan- 
tifying the amount of entanglement that was actually generated can, in principle, be accomplished 
using full-state tomography. This method requires the determination of a parameter set that is 
growing exponentially with the number of qubits and becomes infeasible even for moderate numbers 
of particles. Non-trivial bounds on experimentally prepared entanglement can however be obtained 
from partial information on the density matrix. As introduced in [K.M.R. Audenaert and M.B. 
Plenio, New J. Phys. 8, 266 (2006)], the fundamental question is then formulated as; What is the 
entanglement content of the least entangled quantum state that is compatible with the available 
measurement data? 

O ' We formulate the problem mathematically [J] employing methods from the theory of semi-definite 

programming and then address this problem for the case, where the goal of the experiment is the 
I I , creation of graph states. The observables that we consider are the generators of the stabilizer group, 

thus the number of measurement settings grows only linearly in the number of qubits. We provide 
analytical solutions as well as numerical methods that may be applied directly to experiments, and 
compare the obtained bounds with results from full-state tomography for simulated data. 
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c5 ■ I. INTRODUCTION 

^ . 

/ Detecting and quantifying entanglement is one of the core problems in quantum information theory [3, Q . The 
detection of entanglement can in principle be accomplished by measuring the complete quantum state and, thereafter, 
' applying separability tests. However, the dimension of the density matrix grows exponentially with the number of 
K*" [ constituents of the system. Therefore, full state tomography becomes very costly and experimentally infeasible already 
■ for a moderate number of particles. Thus, it is of interest to detect and quantify entanglement, even when only partial 
' information on the density matrix is known. Entanglement Witnesses represent one way to verify the existence of 
. entanglement with only a few measurements In particular, witnesses for the detection of stabilizer states have 
• ' been constructed in which also provide lower bounds on the fidelity. 
, In it was demonstrated that one may define an entanglement measure on the basis of witness operators and 
Q.^ ' provide lower bounds on entanglement measures 0, Si ■ 

I However, the restriction to witness operators is unnecessary and may neglect information that is obtained when 
i measuring the local operators into which the witness operator has been decomposed A direct calculation of the 
. !^ ' least amount of entanglement (in accordance with an entanglement measure of choice) that is compatible with the 
, measured data of arbitrary observables is proposed in [p and this approach is guaranteed to deliver the best lower 
$_i ' bounds that can be obtained from the information that is available. The same philosophy may also be followed when 
5^ , bounding other quantities such as the fidelity with a desired target state. We apply the method described in [l| to 
the case where the goal of the experiment was the creation of cluster states. The observables we consider are the 
generators of the stabilizer group (for an introduction to the stabilizer formalism see e.g. @). Thus the number of 
measurement settings grows only linearly in the number of qubits. 

This article is structured as follows: First, we describe how to provide lower bounds on the fidelity with a target 
state (here stabilizer states) in Sec. |lT] and in Sec. IHII apply this to some simple examples. Sec. IIVI discusses the 
general approach to estimate robustness measures from incomplete information on the density matrix. Then, in Sec. 
|V|we utilize this approach to obtain lower bounds on the Global Robustness of Entanglement, and give closed formulae 
for systems consisting of two, three, and four qubits. A comparison of the obtained bounds with exact values for noisy 
cluster states is provided in Sec. IVll 
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II. MINIMAL FIDELITY AND ENTANGLEMENT 



111 many experiments we aim at creating a particular pure quantum state |(/)). Needless to say, experimental 
imperfections and noise will usually lead to a noisy approximation to this state, i.e. a fidelity that is different from 
unity. This naturally raises the question as to how close we actually are to the target state. It is desirable to find 
simple sets of measurements that give us enough information to find useful lower bounds on the fidelity that has been 
achieved in the experiment. This problem may be solved with the methods that have been developed earlier in [II]. 
More formally, we will measure a set of obscrvablcs {A^} and find measured mean values a^. Then we will find the 
state p that predicts the mean values Ui and that has the least fidelity with the target state \(f>). Mathematically this 
is formulated as 

Fmtn = min[tr[\(j)){(j>\p] : tr[Aip] = ai,p>0]. (1) 

The solution to this problem is called the primal optimal. This problem is in fact numerically very efficiently solvable 
as it is linear program which is a special case of a semi-definite program. As such there are firstly very efficient 
numerical algorithms and, employing the concept of duality, one can also find lower bounds on the minimization 
problem [To| . Indeed, by duality we find that 

F,„„ = TOax[^A,a,; : -^A,A,) > 0]. (2) 

i i 

The solution to the latter, dual, problem is the dual optimal. That we really have equality, as we have implied 
here, is not trivial, but is true for linear programs, and for general semi-defiiiitc programs it is true under very mild 
conditions (see for details) and is usually safe to assume at the beginning (though that this needs to be checked 
should be remembered when primal and dual optimum do not appear to coincide). The same methods can also be 
used to verify and quantify entanglement measures. In the bi-partite setting this can for example be done for the 
logarithmic negativity 0, [HI, and simple analytical formulae can be given [I| for useful sets of observables. 

For multi-partite settings we have a variety of measures available [p, Q . One may chose simple generalizations of 
the negativity measures but may also study robustness measures [14| . The method for the description of negativity 
measures can be deduced directly from In this note we will thus only present the basic approach for robustness 
measures. 



III. EXAMPLES: FIDELITY ESTIMATION 



A. GHZ States 



The optimization problem formulated above may look somewhat daunting. Let us therefore consider some examples. 
First, we consider the quantitative verification of the fidelity with the EPR state \(f>2) = "^(|00) + |11)). Let us measure 
the expectation values of the observables Ai = Xi (S) X2 and A2 = Zi (g) where X (Z) is the Pauli x (z) operator. 
The unit trace condition on the density matrix is tr[p\ so that we have A3 = t. Then we find 

_ 0^ + 02 

^ mm — ^ • \^ J 



This is seen from the choice 



1 + 02 2ai + a2 - IN 

= 1 \ ^ I - 02 1 - a2 I 

^ ^ 4 I I - a2 I - 02 I ^ ' 

y2ai +02-1 I + a2 



and for the dual problem with the choice Ai = A2 = | and A3 = 0. The verification of the GHZ-fidelity, that 
is the overlap with the state \<j>) = i(|000) + |fff)) may also be considered. Here wc measure the observables 
Ai = Xi ^ X2 X3, A2 = Zi iSi Z2 iSi t, and ^3 = 1 (g) Z2 ® Z3. Again, A4 = 1. Here it is a little harder to find the 
closed formula but it is actually 



(5) 
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The dual optimal is then of course Ai = A2 = A3 = i and A4 = — i. The optimal p for the primal problem has 
Pi,i = max[{a2 + O3)/4,0], pi^j = pi^i = pj,j ~ pj^i for i,j e [2, 7] and pi^g such that the above optimal emerges. 

For a general n-particle GHZ-state we measure for example Ai — Xi (g) ■ ■ ■ X„, and Am — Zm-iZm for m > 2, 
to find 

Fmin = max[ °^ ^ °" — ^ , 0] (6) 
It is straightforward to read off the form of the dual optimal from this expression. 

B. Cluster States 

We find the same bounds on the fidelity for cluster states. As a matter of fact, these fidelity estimates are true for 
any observables which generate a stabilizer group. A proof is given in Appendix [Bl 

IV. ENTANGLEMENT MEASURES 

In this section we will discuss the estimation of entanglement measures from tomographically incomplete measure- 
ments . As mentioned, [i| discusses already the logarithmic negativity measures. For example, when one measures 
X ® X and Z ® Z and finds and then as demonstrated in [l[: Emm = raaxi^ ,log{\ax\ + |az|))- and if one 
additionally measures Y ®Y and finds a.y, then Emm = rfiax(^^ log{l + \ax\ + \o.y\ + \az\))- Here we will present the 
approach for the Global Robustness of Entanglement. For bi-partite systems, this is defined as 

E{p) = min[tr[a] : a > 0, p^ + > 0] (7) 

where a must be Hcrmitian and positivc-scmidefinite, and F denotes partial transposition. For many particles, say n, 
a natural extension is 

E{p) = min[tr[a] : cr > O,^'^" + cr^° > 0,Va G {1, . . . (8) 

Again, given some expectation values ir[p74i] — ai, we would then determine 

Emi7i = min[E{p) : tr[pAi) = a^, p > 0]. (9) 

This is again a semi-definite program and is thus rapidly solvable using numerical programs. For analytical work it 
will again be interesting to derive the dual which will allow us to find lower bounds on the above minimization. This 
derivation can be done in the following steps: 

Emin = rain[tr[a\ : tr{pAi) ~ ai, p^" + > 0, p > 0, c > 0] 
= max[min{tr[a{l — ^^'7q°)] — t'l^blg" p] 

a a 

- ^ pi(tr[pAi] - a^) : p, o- > 0} : ?;„ > 0, p] 

i 

= max[min{-^tr[?£'' p] - ^p^(tr[pAi] - a^) : p > 0} : > 0, 1 > ^?7Q°,/i] 

a i a 

= max[Y^ piUi : 1 > ^ ?7a° , ^ ??a° + ^ PiAi < 0, 77a > 0, p] (10) 

i a a. i 

V. APPLICATION TO STABILIZER STATES 

In this section we will utilize the approach described in the previous section to explicitly calculate lower bounds on 
the Global Robustness of Entanglement. Wc will see that the proper choice of observables transforms the optimization 
problem into a linear program, which may be solved analytically as well as numerically using well-known algorithms 
like the Simplex method. 

We assume now that the goal of the experiment was to create cither a cluster state with the associated adjacency 
matrix Fyi or a GHZ state. Then a natural choice for the observables Ai would be the generators Ki of the abclian sta- 
bilizer group. For cluster states the stabilizers (or correlation operators) are given by Ki = XiZj^. where the subscript 
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Ni is to be understood as applying Z-operators to all neighbors of the i-th qubit in the lattice defined by Va- For 
iV-qubit GHZ states the generators of the stabilizer group are: Ki = Xi® ...®Xi^ and for k = 2, N: = Zk_iZk. 
Due to the commutation relations fulfilled by these operators, it is easy to see that the symmetries that leave the 
observables = tr{pKi) invariant are given by the transformation p — > p' = jtv i^^o -^i^ ■■■^n' P^i^ ■■■^n' ■ 

We may therefore restrict our attention to states of the form: 



P= E c,,...,.Ki-K^ (11) 

ii,...,iN—0 



with real coefRcients Ci^ .. ^^. We can further restrict the matrices rja of the dual problem to have the same sym- 
metries as the states p. The eigenvalues of p are: Aji...j„(p) = J^l^ jjv=o(~l)'^''^ ■■■(~l)*""'"'-»i - »iv > where the 
jk G [0, 1] form a binary index for A. The symmetries obeyed by p also imply that the (unnormalized) state 
a has the same symmetries as This can be seen as follows: One may define the completely positive map 

M"^) = W'J2ii {^=0^1^ ■■■-^n'^-'^i^ ■■■-^n' ■ Then, assume we found i?mi„ and the corresponding operator a, 
such that {p + (t)^° > Va G {1,...,N}. Since A((p + ct)^°) = p^" + A((t)'"°, one concludes that a must 
be invariant under rotations of the stabilizer group. Thus: a — J2i i =o^ii--iN^i^ ■■■^n' j with eigenvalues 

A. 2 Qubits 

1. 2 Qubits: Upper Bound 

In this section we will evaluate the entanglement for a composite system of two qubits, supposedly prepared as a 
cluster state. The measurements performed on this system result in ai = tr[pKi) and 02 = tr(pK2), with Ki = X®Z 
and K2 = Z ® X. W.l.o.g. we restrict to the case of positive ai but write the solutions of the more general case of 
arbitrary a^. The primal problem reads 

min{tT{a) : {p + af' > 0, p > 0, cr > 0, tr{pK,) = a,] (12) 

We denote the eigenvalues of p resp. of its transpose by: 



E i~iy''H-iy'''c,,,, (13) 

ii ,12—0 
1 

A,u.(/0= E (-l)''''(-l)''''("l)^'^'C'i'- (14) 

^1,12— 

Eigenvalues of cr and are of the same form, and we denote the corresponding coefficients with di-^i^ ■ The coefficients 
Coo = 1/4, cio = ai/4, coi =02/4 are given by normalization and measurement constraints. In the case ai + 02 < 1, 
one may set cn = 0, which is the coefficient that changes sign under partial transposition. Thus, cr = in this case. 
Otherwise, upper bounds on tr(a) = 4doo are obtained by the choice doo = (o-i + 0,2 — l)/4, di-^i^ = — (ioo/3 else, and 
cii = — c^oo- The upper bound on the Global Robustness of Entanglement is thus given by: 

Emin = max{0,\ai\ + \a2\ ~ I}- (15) 
2. 2 Qubits: Lower Bound 

Here we will derive a lower bound on the Global Robustness of Entanglement according to Eq. [TOl We will see that 
this lower bound coincides with the upper bound derived in the previous section. First, one may restrict the matrices 
rja to have the same symmetries as p: 

Vc.= ^Ik^'K^- (16) 

11,42=0 
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Since the partial transposes Fi and r2 have the same impact on the rja (both change the sign of the coefficient cn), 
we may simphfy the problem by setting rji = r]2 = 77/2. Again we consider only the case of positive a^'s. The dual 
problem can now be formulated as the following eigenvalue problem in the form of a linear program: 

max[fio + ^ifli + ^202 : yUol + l^iKi + IJ.2K2 + rf' < 0, 77 > 0, i]^^ < 1]. (17) 

Besides the trivial solution (all variables equal zero), a little thought shows that the above system of inequalities is 
fulfilled by — /io = /ii = /i2 = 1 and coo = — cio = — cqi = cn = 1/2. Thus: 

Emin = max{0,\ai\ + \a2\ - I}. (18) 

which coincides with the upper bound presented in the previous subsection. 

B. 3 Qubits 

If the goal of the experiment was the creation of a triangle cluster state, the observables are naturally Ki = X^Z(^Z, 
K2 = Z ® X ® Z ^ and K-^ = Z ® Z ® X with measurement outcomes ai = tr{Kip), i £ {1, 2, 3}. Then, one finds a 
solution similar to the 2-qubit case, in the sense that it only depends on the two largest measurement outcomes: 

Emin = maxiO, \ai\ + \a2\ + \as.\ - min{\ai\, \a2\, \as.\) - 1}. (19) 

C. 4 Qubits 

Let us now consider the case, where the goal of the experiment was the creation of a 4-qubit cluster state associated 
with a square lattice graph (box cluster state). Then the four generators of the corresponding stabilizer group are 
given hyKi=X(g)Z(g)l(g)Z,K2 = Z(^X(g)Z(S)l, K3 = t®Z®X(g)Z,K4 = Z®l®Z®X. The measurement 
outcomes are denoted by = tr{pKi) for is {1, ... ,4} and are assumed to be non-negative. Thus, the problem for 
the square lattice case reads: 

4 3 

max[po + ^ pitti : ^ + ^ fiiK^ < 0, 77^ > 0, ^ 7/^° < 1] (20) 

i=l i=0 

where Kq ~ 1. The matrices 7]a are restricted to: rja — X^ii ii=i)^i^ ij^i^ ■ ■ -^l"- The partial transposes are 

therefore given by: 77^= = En,...,.,=o(-l)*^ '''-^...J<^ ■ • ■ 

Even though, this translates to a system of inequalities which looks rather complex, one may realize easily that 
the solution of the two-qubit case represents also a solution for this system. This means, /io = — 1, Mi = M2 = 1 and 
A*3 = A*4 = 0, and regarding the 77^ one obtains Cqooo = c^\|,(, = \, cJo|,o = ^ -\ and 7/2 = 77i> % = = 0. One 
may easily check, that another solution is given by the following set of parameters: /^q = ~5, /ii =2, /i2 = 2, /i3 = 
2, /i4 = 2. and the coefficients of the operator sum representation of the r\a are listed in Tab. [Ill Summarizing these 
results gives: 

Era%n = max^, (|ai| + \a2\ - 1), 2(|ai| + |a2| + \a-i\ + |a4|) - 5] (21) 

VI. QUALITY OF THE ESTIMATE AND LOCAL STATISTICS 

In order to check the usefulness of the obtained bounds, we compare these bounds with exact values for simulated 
noisy cluster states. We assume that after a perfect cluser state was created, the qubits are subject to local dephasing 
for a certain time (here we assume 10 ms). Then, the system is described by the following master equation: 

P=\Y.{Z,,pZ,,-p) (22) 

i 

where 7 is the dephasing-ratc. which we take to be (10s)~^. A comparison between exact values of the Global 
Robustness with our estimate is given in Tab. [J It shows that the estimate deviates only a few per cent from the 
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No. Qubits 


exact value 


estimated value 


relative deviation 


2 


0.8142 


0.8097 


0.0055 


3 


0.8185 


0.8097 


0.0108 


4 


2.2995 


2.2387 


0.0264 



TABLE I: Comparison of the Global Robustness of Entanglement (GRE) for 2, 3, and 4 qubit noisy cluster states with an 
estimate of the GRE from measurements of the generators of the stabilizer group only 

exact value. It is obvious that the bounds can be improved by considering any additional information on the density 
matrix. When one performs measurements on distant parties, the observables such as {X ® X) must be gained from 
local measurements. The entanglement depends mainly on the correlations, as we have seen in the above comparison. 
However, the obtained local statistics may be used to improve the bounds. Consider the measurements {X®X) = 0.9, 
{Z®Z) = 0.7, {Z®t) = 0, (1 = 0.25. In this example, the GRE yields 0.6 when one considers only the XX- and 
ZZ-observables. Taking into account the local statistics improves the GRE by more than ten per cent to 0.6671. This 
example shows that one should use all available information the measurements provide to obtain optimal bounds on 
entanglement. 

VII. CONCLUSION AND OUTLOOK 

In conclusion, we investigated how to obtain lower bounds on the fidelity and robustness measures from partial 
information on the density matrix of a multi-partite system. We utilized the symmetries of the stabilizer group to 
formulate the problem as a linear program, which can be treated analytically as well as numerically. Analytical 
solutions were obtained for two, three, and four qubit-sy stems. This method is of particular interest for experiments, 
since the number of measurement settings grows only linearly in the number of qubits, whereas full-state tomography 
requires an exponential number of settings. A comparison of the obtained bounds with exact values of the Global 
Robustness shows that the difference is in the order of only a few per cent. 

For the future, it will be interesting to investigate if analytical solutions to our approach can be found for systems 
with an arbitrary number of constituents. 
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APPENDIX A: SOLUTION TO THE DUAL PROBLEM FOR A 4 QUBIT MEASUREMENT 



a 


1 


2 


3 


4 




1 


2 


3 


4 


J") 

^0000 


3/16 


3/16 


3/16 


3/16 


''0110 


-1/16 


1/16 


1/16 


-1/16 


(a) 
'^1000 


-1/16 


-2/16 


1/16 


-2/16 


''0101 


1/16 


-3/16 


1/16 


-3/16 


(a) 
•^0100 


-2/16 


-1/16 


-2/16 


1/16 


''0011 


-1/16 


-1/16 


1/16 


1/16 


(a) 
"^0010 


1/16 


-2/16 


-1/16 


-2/16 


''1110 


2/16 


-1/16 


2/16 


1/16 


(a) 
^0001 


-2/16 


1/16 


-2/16 


-1/16 


''1101 


-1/16 


2/16 


1/16 


2/16 


^1100 


1/16 


1/16 


-1/16 


-1/16 


'-1011 


2/16 


1/16 


2/16 


-1/16 


J") 

■-lOlO 


-3/16 


1/16 


-3/16 


1/16 


''0111 


1/16 


2/16 


-1/16 


2/16 


J") 

^1001 


1/16 


-1/16 


-1/16 


1/16 


^(«) 

''1111 


-1/16 


-1/16 


-1/16 


-1/16 



TABLE IL Coefficients for the operator sum representation of the operators r]a, where the goal of the experiment is a box 
cluster state 



APPENDIX B: PROOF FOR THE FIDELITY OF NOISY CLUSTER STATES 

Let ai ~ tr{Kip), i £ {1,...,N} be the mean values of the stabilizer operators. Twirling over 
the stabilizer group allows us to restrict to states p — ^^J" ■■■ ijv=o ''^i - 'w-^'^'r ' ' ' ^./v*- with eigenvalues 
•^ji----,ijv (p) = 2^ X^ii •■■ ijv=o(^-'^)^"' Furthermore, the target state may be written as: = 

Primal problem: Now we choose the coefficients Ci-^...iN = Sa^i ~ X)*fc + 1- One may assume, that the mea- 
surement outcomes Ui are positiv. Any other choice could be achieved by local rotations of the qubits. Then, p > 
follows from the fact that p is diagonal in the stabilizer basis, and all coefficients are positive. Because V m : K^^^ = 1, 
we find 

tri\^){<f>\p) = ^ E ^'1--" (Bl) 

ii,...,iN—0 

IN 

= ^ E (E*'»«^-E*^ + l)=2(E«fc-^ + 2) (B2) 

2i,...,?jv— k—1 k 

Dual problem: The dual problem may be solved by the choice Aq = —N/2 + 1 and A.; = 1/2 for z > 1. In order to 
check the validity of this solution, one must prove that x — S > with x — \4'){4'\ a-^"^ = h J2iLi Ki + (1 — y)-"^- 
The eigenvalues of S are given by 



A,....,«(S) = -f + l + (B3) 



1=1 



This means Ao...o(S) = 1. Since Ao...o(x) = 1; we have Ao...o(x — S) = 0. It is easy to see that Aji...j„(^) < for 
(ji, . . . , jjv) 7^ (0, ... , 0), and Aji...j„(x) = for these indices, thus x — 2 > 0. 



